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In this talk, I describe the general characteristics of radiative neutrino mass models that 
can be probed at the LHC. I then cover the specific constraints on a new, explicit model 
of this type. 
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1. See-saw versus radiative neutrino mass generation 


Many models of radiative neutrino mass are possible, some of which have been anal¬ 
ysed in depth, while others have been examined only briefly, with quite a number 
yet to be explicitly written down let alone analysed. A useful organising principle 
for Majorana neutrino mass models is to use standard model (SM) effective opera¬ 
tors that violate lepton number conservation by two units {AL = 2 operators) as a 
base from which to construct new theories.— “— 

Restricting ourselves to operators containing SM fermions and a single Higgs 
doublet, they occur at odd mass dimensions, with the lowest order one occurring 
at d = 5. That operator is the famous Weinberg operator, which has the schematic 
structure LLHH, where L is a lepton doublet and H the Higgs doublet. After 
electroweak symmetry breaking, this operator (actually a set of operators because 
of the family structure) directly induces a Majorana neutrino mass for left-handed 
(LH) neutrinos given by the see-saw formula. 



( 1 ) 


where M is the scale of the new physics that gives rise to the operator in the 
low-energy limit. 

Underlying renormalisable theories yielding LLHH are constructed by “opening 
up” the operator. The type-1,-2 and -3 see-saw models are the minimal, tree-level 
ways to open up LLHH By starting with the effective operator, one may 
systematically construct all of the minimal underlying models. 

Other AL = 2 effective operators have mass dimension seven or higher and 
feature other lepton and quark fields in their expression, except for those in the 
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generalised Weinberg class LLHH{HH )'^These additional particles have to 


be closed off in loops to produce a Majorana neutrino mass self-energy graph from 
the operator. Thus, theories based on c? > 7 operators (apart from the generalised 
Weinberg class) necessarily produce radiative neutrino mass generation. Underlying 
renormalisable theories can then be systematically constructed by opening up the 
operators in all possible ways, subject in practice to minimality assumptions, just 
as the three see-saw models may be derived from the Weinberg operator. 

A list of all gauge invariant AL = 2 operators at c? = 5, 7,9,11, constructed from 
SM fermions and one Higgs doublet, was produced by Babu and Leung (BL).— The 
ninth operator in their list, Og = LLLe'^Le'^ (which has d = 9), is the basis of 
the historically important Zee-Babu model of neutrino mass generation. Two 
exotic scalars are introduced: an isosinglet, singly-charged state h coupling to LL 
and a doubly-charged, isosinglet k coupling to Both carry two units of lepton 
number. The cubic term hhk combines with those interactions to induce AL = 2 
and produce Majorana neutrino mass at 2-loop order. ATLAS has searched for k 
through the same-sign dilepton channels ee, e/i and /r/i, deriving a lower bound of 
about 320 GeV on the mass^ (see also^). 

It is nice that the AL = 2 operator perspective places the radiative neutrino mass 
models at one end of a systematic list of Majorana mass models, bookended by the 
tree-level see-saw models at d = 5. From the phenomenological viewpoint, radiative 
models are interesting because the mass scale of the new physics is lower than the 
favoured very high scale for see-saw models, and has more chance of encroaching 
into the LHC regime. 

It is very much worth noting that upper bounds on the see-saw scales can be 
derived from naturalness considerations, because the new particles and interactions 
will destabilise the electroweak scale if the new physics occurs at too high a scale. 
The type-1 upper bound was first computed by Vissani— to be 3 x 10^ GeV, if 
the Higgs self-energy graph from neutrino Yukawa interactions is not to produce a 
correction to the parameter in the Higgs potential greater than (1 TeV)^. This 
one-family result has recently been generalised to the realistic three-family case,— 
producing the bounds 



Mn, < 4 X 10^ GeV, Mn^ < 7 x 10^ GeV, MjVj < 3 x 10^ GeV 


where Ni^ 2,3 are the three heavy neutral lepton mass eigenstates with < Mn^ < 
MjVg, and mmin is the minimum light neutrino mass. It is interesting that these 
bounds imply that standard thermal, hierarchical leptogenesis^, which requires 
sufficiently massive heavy neutral leptons, must involve some level of fine-tuning for 
the minimal type-1 see-saw model. This is true for all three cases: Ni, N 2 and N 3 
leptogenesis. Naturalness, if successful leptogenesis is desired, may be restored by 
extending the theory. An obvious example is the supersymmetric extension, and 
another is to add extra Higgs doublets, to divorce the leptogenesis parameter space 
from neutrino mass generation (one version of this is discussed in Ref.—). 
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2. Opening up d — 7 operators 

Let us examine the AL = 2, d = 7 operators from the BL list as the next most 
complicated cases after the Weinberg operator. The flavour contents are 

O 2 = LLLe'^H, O 3 = LLQd^H, O 4 = LLQu'^H, Og = Le‘'u^d^H, (3) 

where the BL numbering scheme has been used, with Q the quark doublet, d the 
RH down quark and u the RH up quark. Flavour structures 03,4 each yield two 
independent operators once weak-isospin index contraction possibilities are taken 
into account. Operator O 3 is the basis of the pioneering 1-loop Zee model,— while 
models constructed from O 3 and Og have been analysed in depth by Babu and 
Julio.— In addition, d = 7 contains the Weinberg-operator generalisation O 4 = 
LLHHHH.^ 

We adopt the minimality assumption that the underlying renormalisable theories 
obtained from opening up these operators contain only exotic heavy scalars and 
vector-like fermions. Table [T] lists the quantum numbers under SU(3)xSU(2)xU(l) 
for scalar-only and scalar-fermion extensions.— The example in boldface will be 
analysed in detail below. The diagram topologies are given in Fig. [TJ 


Table 1. Quantum numbers of exotic scalars and 
fermions in underlying theories for the d = 7 oper¬ 
ators. 


Scalar 

Scalar 

Dirac fermion 

Operator 

(1,2,1/2) 

(1,1,1) 


02,3,4 

(3,2,1/6) 

(3,1,-1/3) 


03,8 

(3,2,1/6) 

(3,3,-1/3) 


03 

(1,1,1) 


(1,2,-3/2) 

02 

(1,1,1) 


(3,2,-5/6) 

03 

(1,1,1) 


(3,1,2/3) 

03 

(3,2,1/6) 


(3,2,-5/6) 

03 

(3,1,-1/3) 


(3,2,-5/6) 

03,8 

(3,3,-1/3) 


(3,2,-5/6) 

03 

(3,2,1^) 


(3,3,2/3) 

03 

(1,1,1) 


(3,2,7/6) 

04 

(1,1,1) 


(3,l,-l/3) 

04 

(3,2,1/6) 


(3,2,7/6) 

Os 

(3,2,1/6) 


(l,2,-l/2) 

Os 

(1,4,3/2) 


(1,3,-1) 

o[ 


It is clear that searches for radiative neutrino mass models are included in gen¬ 
eral searches for exotic scalars and fermions. Each specihc model has a particular 
scalar-scalar or scalar-fermion pair. In addition, the decay branching ratios for the 
exotics are typically quite constrained by the requirement to reproduce the correct 
neutrino masses and mixing angles. This additional information must be used when 
interpreting the experimental constraints, and it constitutes a “smoking gun” for 
exotics that are responsible for radiative neutrino mass generation. 
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(a) (b) (c) 


Fig. 1. Diagram topologies for opening up d = 7 operators, (a) Scalar-only, (b) Scalar plus 
fermion for 02,3,4,8- (c) Scalar plus fermion for O'^. 

3. Collider constraints on a new radiative model 

We now look at the boldfaced case in Table [T]— Our theory will produce O 3 , with 
a subdominant Og. The exotic scalar ^ and exotic fermion x have the quantum 
numbers 

0^(3,1,1/3), (3,2,-5/6). (4) 

The Lagrangian is 

- £ = + h.c.) 

+ (Y^f^L^Q.cj, + + Yf^^d,xH + h.c.) (5) 

where we use two-component notation for the fermions and i,j are family indices. 
The Yukawa coupling constants play no role in neutrino mass generation and 
are thus set to zero for simplicity. We also impose baryon-number conservation to 
forbid the QQ4>^ and ducj) interactions permitted by the gauge symmetry. 

The diagram responsible for neutrino mass generation is given in Fig. [21 It is 
proportional to the down-type quark masses. Barring anaesthetic hierarchies in 
the coupling constants, the terms proportional to the 6 -quark mass will dominate, 
which is what we will assume is the case. Associated with that, for simplicity we 
switch off the mixing of x with first- and second-generation quarks. 


r 1 Qi d X X \ T 

Ln —►-*—« y >■—f ^ X >•—*-<— Lif 


{Hf \h) 


Fig. 2. Feynman diagram responsible for 1-loop generation of neutrino masses and mixings. 


The neutrino mass matrix is then 




3 

IGtt^ 



+ (* ^ j)) 


rribmB , rn% 

mbB^2 - 2 " — T- 

mi - ml 


( 6 ) 
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Fig. 3. Lepton flavour violation bounds, (a) Allowed region in the |(^|, plane for m-^ = 2 TeV. 
(b) Allowed region in the |^|, plane for = 2 TeV. 


The bottom quark mixes with the isospin +1/2 component of x to form two mass 
eigenstates, b and B, with masses mb and ms respectively. The parameter mbs = 
V/\/2^ where {H^) = v/\/2. At this level of approximation, there are two 
massive neutrinos and a massless one. In reality, the lightest neutrino will pick up a 
small mass, but its value is unimportant and will be neglected. CP violating phases 
will also be set to zero for simplicity. 

The rank-2 mass matrix may be expressed as the symmetrised outer product of 
two vectors a±, 

mi, = a+OL^ + a-O^i (7) 


where 


NO 


>±i 


= {y/^U2 ± iy/^ul) 


af 


± iy/m2U2) , 


( 8 ) 


with NO and 10 denoting normal and inverted ordering, respectively. The vec¬ 
tors ui. 2,3 are the columns of the PMNS matrix, and the complex number is a 
Casas-Ibarra-like parameter, not determined by low-energy experiments. Equating 
m and [7] constrains the parameter space to the region compatible with the neutrino 
oscillation results. 

The next set of constraints on the parameter space come from the lepton flavour 
violating processes /r —>■ ey, /i —>■ eee and /iiV —>■ eN An example of the results is 
presented in Fig.|3l The green dash-dotted, blue dotted and red-dashed lines bound 
the regions corresponding to BR(/r —ey) < 5.7 x 10“^^, BR(/r —>■ eee) < 10“^^ 
and BR(/rAu —eAu) < 7 x 10“^^, respectively. The grey region is thus excluded. 
The magenta dashed line shows the reach of the /rTi —>■ eTi Mu2E and COMET 
experiments. The allowed region is divided into B and T regions in which BR(^ —>■ 
hv) = 1 and BR((/) —>• 6z/) < 1, respectively. 

The vector-like fermion doublet x consists of the charge —1/3 quark B' which 
mixes with the bottom quark, and an exotic charge —4/3 quark Y. CMS have 
searched for R-like particles, whose dominant decay modes are B ^ Zb and B Hb 
in our model with branching ratios that are strongly constrained by the need to ht 























April 10, 2015 0:6 


WSPC Proceedings - 9.75in x 6.Sin 


Volkas page 6 


6 


the neutrino oscillation data. Taking that into account, the CMS bound is ms > 620 
GeV.— No Y search has been performed. 

The other collider constraint comes from searches for leptoquark scalar (j). It is 
pair-produced at the LHC from gluon-gluon fusion and qq annihilation. Since the 
production is through its colour charge, the rate depends only on m^. For example, 
for = 500(600) GeV, the production cross-section a{pp is 82(23.5) fb. 

The main decay modes are into Lt and biy, where L = (e, p, r). In the parameter 
region my^B ^ so that LY and Bv final states are kinematically forbidden, 
the partial decay rates are. 


T{(j) Lt) 
r((/) VLb) 


rritj, 

Stt 

Stt 


Y, 


LQ4, 


L3 


Y^^^ COS 02 


2 


r, 


Lx4> 


sine/i 




(9) 


where / is given in Eqs. 4.16 and 4.21 of Ref.—, and the mixing angles 0 i^2 are 
defined through 


sin 01 


mbBm^ 

9 2 ’ 

— mp 


sin 02 


rribBrnb' 

9 2 ’ 

— mp 


( 10 ) 


where my = ybv/\/2 is the bottom quark mass in the absence of mixing with B'. 
The decay rates depend on the same Yukawa coupling constants that also contribute 
to neutrino mass generation, so are quite constrained, as well as ICj. 

Recall that the allowed regions B and T in Fig. [3] correspond to regimes where 
BR((/) —>■ bi/) = 1 and BR((^ —>■ bv) < 1, respectively. In the B region, the main 
signature is thus pp —>• (j)(t> bb + missing Ey , where sbottom pair searches apply. 
The bound in this case is nirf, > 730 GeV at 95% C.L.— The analysis for region T 
requires recasting stop searches as well and is summarised in Fig.|Tl which plots the 
various branching ratios as functions of m^, as well as the limits from ATLAS and 
GMS.— In numbers: The lower limit from GMS in the 66-|-missing Ey channel 
is in the range 520 — 600 GeV. From the (e, p) + missing Ey + (6—)jets, ATLAS sets 
a limit of approximately 580 GeV. Finally, the {e, p)~^ {e, p)~ -I- missing Ey -|-jets 
channel yields about a 620 GeV limit from ATLAS data. 


4. Final remarks 

A systematic procedure for constructing radiative Majorana neutrino mass models 
has been developed from standard model AL = 2 effective operators containing 
leptons, quarks and the Higgs doublet. The LHG is a useful tool for searching for 
the exotic scalars and fermions that appear in these models, and specific bounds for 
a new model were summarised in this talk. The 13 TeV LHG will extend the reach 
of past searches. Some naturalness concerns for hierarchical, thermal leptogenesis 
in the minimal type-1 see-saw model were also mentioned in passing. 
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(GeV) 


Fig. 4. Decay branching ratios of 0 as a function of rrifj) in region T. The curves with positive 
slope running to the top right show the ATLAS (thin line) and CMS (thick band) limits. 
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